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Abstract 

Residence times of a particle in both the wehs of a double-well system, under 
the action of zero-mean Gaussian white noise and zero-averaged but tempo- 
rally asymmetric periodic forcings, are recorded in a numerical simulation. 
The difference between the relative mean residence times in the two wells 
shows monotonic variation as a function of asymmetry in the periodic forc- 
ing and for a given asymmetry the difference becomes largest at an optimum 
value of the noise strength. Moreover, the passages from one well to the other 
become less synchronous at small noise strength as the asymmetry parameter 
(defined below) differs from zero, but at relatively larger noise strengths the 
passages become more synchronous with asymmetry in the field sweep. We 
propose that asymmetric periodic forcing (with zero mean) could provide a 
simple but sensible physical model for unidirectional motion in a symmetric 

periodic system aided by a symmetric Gaussian white noise. 
PACS numbers: 82.20.Mj, 05.40.+j, 75.60.Ej 
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Several physical models have recently been proposed |T|-0,3 to understand possible av- 
erage asymmetric motion of a Brownian particle in a periodic potential. Living systems 
are manifestly nonequilibrium and quite understandably such an asymmetric motion has 



work out of nonequilibrium systems is not new, the biological experimental observation has 
given enough motivation recently to renew effort in that direction. It has resulted in better 
understanding of the problem and also it has helped in inventing new devices for practical 
use In the present work we study a symmetric two-well system subjected to zero-mean 
Gaussian white noise. We apply an external field that is periodic in time. The external 
field is taken to be temporally asymmetric but with mean force zero in a period. We find 
that eventhough the mean deterministic force experienced by a particle due to the external 
field is zero, the Gaussian white noise (centered at zero) helps it to extract work while mov- 
ing on a potential surface. This is reflected in the asymmetric passages between the two 
symmetrically connected wells modulated periodically by the external field. 

We consider the symmetric two-well potential represented by f/(m) = — + and 
consider the external field h{t) to be periodic in time and assume an asymmetric saw-tooth 
form for it. The asymmetry in the saw-tooth form comes because the positive and negative 
slopes are taken to have different magnitude. The mean force in a period (because of the 
external field) is assured to be zero by taking the maxima and minima of the saw-tooth to 
have values +ho and —ho, respectively. A Brownian particle will thus experience a combined 
(time dependent) potential 



We consider the time evolution of the particle coordinate m{t) to be governed by the over- 
damped Langevin equation 



been observed recently in biological systems ||^. Though the quest for extracting useful 
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where f{t) is a randomly fluctuating force and is taken to be Gaussian with statistics 
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Here < ... > represents average over a large number of realizations of the random forces. 

Our calculation involves solving the Langevin equation numerically and to monitor the 
time evolution of m(t) for a long time for given noise strength D. The calculation is done for 
a fixed subcritical hit) with amplitude ho < he, where he is the minimum value of \hit)\ at 
which one of the two wells of $(m) becomes unstable and disappears. Since hit) < he always, 
the barrier between the two wells never vanishes. Therefore to pass from one well to the 
other a particle need necessarily have to surmount a nonzero potential barrier and therefore 
has to be noise aided. Along the time axis we record the events of passages between the two 
wells. We consider passage to take place from a given well as and when the trajectory m(f), 
emerging from the given well, crosses the inflexion point on the other side of the potential 
barrier separating the two wells. Prom the markers recorded on the time axis we obtain the 
distribution pi(r) of residence times, r, in the well 1 and (similarly for the well 2). And 
also the distribution puih) of field values hit) at which passages take place from the well 
1 to the well 2 (and similarly from the well 2 to the well 1) are calculated from the same 
recordings. 

The distributions puih) and P2iih) determine the evolution of the fraction of population 
in a well as the external field hit) varies. For example, the fraction m2(/i) of the population 
in the well 2 evolves [from (n — l)th step to nth step] as 

m2(n) = m2(n - 1) - m2(n - l)p2lihn-l)ihn - hn-l) + mi(n - l)pi2(/in-l)(^n - ^n-l), 



where /i„ is the field value at the nth subdivision point in a cycle of h. The interval of 
uniform subdivisions (/i„ — /in-i) are taken to be optimally small for better accuracy. In 
our calculation we take (/i„ — hn-i) — Ah — 0.001/ic- So the whole period is divided into 
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= equal segments. This evolution equation together with the periodicity condition, for 
instance, m2(n = 0) = m2{n = ^) gives the closed hysteresis loop Th(h) = m2{h) — mi(h). 
Also, throughout our calculation, we take a = 2.0, and b = 1.0. The hysteresis loop area is a 
good measure of degree of synchronization of passages between the two wells. For example, 
if the passages take place only when the potential barrier for passage is the least, i.e., at 
h = ±/io, the distributions puih) and P2i{h) will be sharply peaked at h = and h = —Hq, 
respectively. In this case the hysteresis loop will be nearly rectangular and therefore will 
have the largest area. On the other hand, if the passages take place all over and randomly 
(corresponding to the case of least synchronization) so that puih) and P2i{h) are uniform 
the loop area becomes the least. We explore how the hysteresis loop area A changes as a 
function of the asymmetry of the field h{t) and also as a function of the noise strength D. 
The asymmetry of field sweep. A, is defined as A = — Tjp-? where Ti is the time for the field 
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to change from /iq to —ho and Tq is the period of oscillation of h{t). Fig. 1 shows a typical 
hysteresis loop with A 7^ as compared to one with A = 0. Notice that the hysteresis loops 
do not saturate for the field amplitude ho = .7 he considered here. 

Figure 2 shows the variation of hysteresis loop function of the asymmetry 

parameter A for fixed values of h^, Tq, and various values of D. The area is the largest at 
A = for small D but is the lowest for relatively larger D. In both the cases, however, 
A{—A) = A{A). This result is to be compared with our earlier work on two-well systems 
1^ where first-passage times instead of residence times were calculated. In that work the 
upper half of the hysteresis loop, corresponging to passages from well 2 to well 1, was 
obtained from the first-passage-time distribution p{h) that spread between /iq and — /iq (and 
not over the whole period) and the other half was obtained by symmetry. Consequently 
the hysteresis loops, by construction, saturated to m(/io) = l.Oh^ and fn^—ho) = — 1.0/ic, 
and were symmetric for all ho including ho = -The- The variation of hysteresis loop area, 
however, showed asymmetry and attained a peak (usually not at A = 0) as a function of A. 
In Fig. 3, we show how, in the present work, area changes as a function of D for fixed A, 
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To, and ho- However, before discussing the physical significance of these results we consider 
the mean residence times fi and f2 in the two wells as a function of A. 

From the distributions Pi(t) and P2{t~) of the residence times we calculate the mean 
residence times in each of the two wells. We, then, calculate the fraction of times, /i and 
/2, the particle spends, on the average, in the two wells. The difference, M = f2 — /i, gives 
a quantity analogous to magnetization (normalized) in magnetic systems. In Fig. 4 we plot 
M as a function of A. From the figure it is clear that M(A) = — M(— A) (upto the order 
of accuracy of our numerical calculation) and vary roughly monotonically. This indicates 
that the particles will tend to accumulate in the well 2 if A > and in the well 1 if A < 0. 
This conclusion is plausible because in the situation under consideration the mean residence 
times fi and f2 add up to f = fi + f2 which is larger than Tq, the period of oscilation of h{t). 
This simply indicates that passages do not take place in every cycle of h{t). For A > 0, 
for instance, there will be larger number of passages from well 1 to 2 than from 2 to 1 per 
cycle of h{t), in an ensemble. This results in a net accumulation of particles in the well 
2, asymptotically. This asymmetry in passages is also refiected indirectly in the hysteretic 
property of the system. For A 7^ the hysteresis loops are asymmetric. Also, as A deviates 
from zero the hysteresis loop area decreases for small D (but increases for relatively larger 
D) [Fig. 2], which indicates that now the passages are less(more) synchronous with the input 
signal h{t). This effect indicates as if a net average constant field is applied in a direction 
determined by A. 

The net accumulation M in a two-well system for given A changes with the noise strength 
D. Fig. 5 shows that M increases initially, reaches a maximum, and then decreases gradually 
as D is increased. We, thus, have an optimum value of D at which the accumulation in the 
well 2 (when A > 0) is the largest after a large number of cycles of h{t). All the results 
obtained in the present work are susceptible to experimental verification by the recently 
developed optical interferometric techniques . 

All the calculated results discussed so far are valid for a double-well potential. However, 
it is not difficult to envisage the situation in case of a periodic potential. In a two-well 



potential, for A > 0, as mentioned earlier, the number of passages taking place from well 
1 to well 2 per cycle of field sweep is larger than the number of passages from well 2 to 
well 1. One may justifiably extrapolate this result to state that in a periodic potential (that 
may even be symmetric), in a given number of cycles of h{t) there will be more 1 — > 2 
passages than 2 — > 1 passages, and hence there will be a net current of particles in the right 
direction (1 — > 2) when A > and in the reverse direction when A < 0. This current will 
increase with the magnitude of A. However, for a given A we can find an optimum value 
of the Gaussian white noise strength D at which the current will be maximum. This is an 
important observation because here we have a physical model for unidirectional motion of 
a particle in a nonratchetlike symmetric periodic potential aided by symmetric Gaussian 
white noise (fluctuating forces). 
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FIG. 1. Hysteresis loops m{h), for = O.lhc, D = 0.15, and Tq = 28.0, axe plotted for (a) 
A = (solid line) and (b) A = 0.5 (dotted line). 

FIG. 2. Plots of hysteresis loop area A versus A for (a) D = O.l(o), (b) D = 0.15(D), and (c) 
D = 0.2(o), (d) D = 0.5(A), and (e) D = 0.7(v)- 

FIG. 3. Plot of hysteresis area as a function of D for (a) A = 0.5(o), and (b) A = j|(D). 

FIG. 4. Shows the variation of accumulation M in a well as a function of A for (a) D = 0.15(n), 
(b) D = 0.2(o), (c) D = 0.5(A), and (d) D = 0.7(v). 



FIG. 5. Plot of M versus D for (a) A = 0.5(o), and (b) A = i|(n). 
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